Oscillations of the Nernst coefficient in bismuth 



O 

(N 

D 
tin 



CD 



o 



00 

m 
O 



X 



Yu. V. Sharlai and G. P. Mikitik 

B. Verkin Institute for Low Temperature Physics & Engineering, 
Ukrainian Academy of Sciences, Kharkov 61103, Ukraine 

We calculate the magnetic-field dependence (oscillations) of the Nernst coefficient in bismuth at 
low temperatures for the case when the magnetic field is directed along the trigonal axis of the 
crystal. In the calculations we take into account the scattering of the electrons and holes in bismuth 
on impurities and the dependence of this scattering on the magnetic field. The results of these 
calculations are compared with the experimental data recently published. 
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I. INTRODUCTION 

In a recent paper, ^ oscillations of the Nernst coeffi- 
cient in bismuth were observed for magnetic fields di- 
rected along the trigonal and bisectrix axes of the crystal. 
These oscillations have the shape of rather sharp peaks 
which originate from the crossing of the Landau levels 
of charge carriers in bismuth with the Fermi level Ep oi 
this semimetal. Moreover, several unusual peaks of this 
coefficient were discovered for very high magnetic fields 
(14 < < 33 T) parallel to the trigonal axis.^ At such 
magnetic fields almost all the Landau levels are empty, 
and one can hardly expect that the unusual peaks re- 
sult from the above-mentioned crossing. So the authors 
of Ref. 2 suggested that the unusual peaks are caused 
by some collective effects in the electron system of bis- 
muth. Interestingly, in the same interval of high mag- 
netic fields several jumps of magnetization were observed, 
which were ascribed to field-induced instabilities of the 
ground state of interacting electrons in bismuth<^ 

Within a simple model of the electron energy spectrum 
of bismuth in magnetic fields, the positions of the usual 
peaks in the magnetic-field axis were explained,—!^ by cal- 
culating the crossing of the appropriate Landau levels 
with Ep. It turned out that most of these peaks are due 
to the hole Fermi surface of bismuth, while some result 
from its electron Fermi surface. However, the unusual 
peaks observed in the Nernst signal at high magnetic 
fields cannot be explained in this way. Nevertheless, we 
showed- that in principle, the positions of some of these 
peaks can be reproduced if one assumes that a small de- 
viation of the magnetic-field direction from the trigonal 
axis occurred in the experiments. In order to support 
or to disprove this possible explanation of the unusual 
peaks, we also suggested^ to investigate experimental in- 
vestigation of the dependence of the Nernst-coefficient 
peaks on the tilt angle of the magnetic field. These ex- 
periments were described in Ref. iGf, and it was found that 
the angular dependences of the usual peaks are in agree- 
ment with the theory^^i^ while the dependences of the 
unusual peaks are essentially different, and hence their 
origin cannot be explained by a small deviation of the 
magnetic-field direction from the trigonal axis. In addi- 
tion, the results of the calculations in Ref. show that 



the usual electron peaks observed in the experiment cor- 
respond to only some of the crossings of the electron Lan- 
dau levels with Ep, while other crossings do not manifest 
themselves. 

In order to continue analysis of the experimental 
findings,— i^i^ in this paper we calculate not only the po- 
sitions of the peaks in the Nernst coefficient but also 
their shape, i.e. the dependence of this coefficient on 
the magnetic field, in the case when the field is along 
the trigonal axis. We carry out the calculations within 
the same model of the spectrum of bismuth as in Ref. i. 
This model is briefly described in Sec. II. In Sec. Ill we 
present formulas for the Nernst coefficient, taking into 
account the effect of impurities on the density of states 
of charge carriers in bismuth, and in Sec. IV we compare 
the results of our calculations with the experimental data. 



II. ENERGY SPECTRUM OF BISMUTH 

The Fermi surface of bismuth consists of one hole ellip- 
soid located at the point T of its Brillouin zone and three 
electron ellipsoids centered at the points L and elongated 
along the bisectrix directions.^ As in Ref. J, we shall de- 
scribe the electron spectrum of bismuth in the magnetic 
field H by the model of Smith, Baraff and Rowell.® Ac- 
cording to this model, the nth Landau level e„ for an 
electron with the quasimomentum pn along H can be 
found from the equation: 



Eil 



p^±\gP,H, (1) 
2mH 2 



where the signs ± correspond to the electron spins that 
are antiparallel and parallel to H, respectively; the en- 
ergy E is measured from the edge of the conduction band; 
LJc is the cyclotron frequency 



\e\H 



Eg is the gap between the conduction and valence bands 
at the point L; g is the effective electron g factor at this 
point; /3o is the Bohr magneton; c is the speed of light; 
e is the electron charge, and the longitudinal and the 
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TABLE I. Parameters of the Smith-Baraff-Rowell spectrum^ 



Electrons 


mil 


m22 




"123 


Orbital mass 


00113 


0.26 


00443 


-0.0195 


Spin mass 


0.00101 


2.12 


0.0109 


-0.13 


Holes 




Ml = M2 




Ma 


Orbital mass 




0.07 




0.69 


Spin mass 




0.033 




200 


Eg = 15.3 meV 






Eo = 38.5 meV 





cyclotron masses, mn and nic, are given by 
rriH — h • m'^ • h, 
nic = [detm^/TOif]^/^. 



(2) 
(3) 

Here h is the unit vector in the direction of the magnetic 
field H. The efi^ective mass tensor has the form: 



mil 
m22 m23 
77123 ni33 



(4) 



where the axes 1 and 3 coincide with the binary and the 
trigonal axes, respectively, while the axis 2 is along the 
bisectrix direction. The effective g factor. 



det 



(5) 



is defined in terms of a spin-mass tensor that has 
the form similar to Eq. Within the model of Smith, 
Baraff, and Rowell^ the elements of nig differ from the 
elements of m'^. 

Since the spectrum of the holes at the point T is 
parabolic, their Landau levels can be easily found^ 



Eq - E = {n+ ^)hujc 



± -gp^H, (6) 
2mH t' ' ^ ' 



where Eq is the edge of the hole band at this point of the 
Brilloun zone. The cyclotron frequency Wc, the masses 
rric and mjj, and the g factor are defined by the same 
formulas ([2]), (jS]), © as for the electrons, but now the 
tensor of the effective masses for the holes m'' has the 
form: 



m 



/ Ml 
A/2 
y M3 



(7) 



and a similar expression is valid for the spin-mass tensor, 
rtij. All the parameters in Eqs. ([H - (O are known for 
bismuthf^i^ see Table HI 

The above formulas lead to the following expression 
for the contribution of the i-th electron ellipsoid to the 
density of electrons in bismuth;^ 



Nt{EF,H) = 



\/E*p - En,; 



(8) 



where 

E*p = EFil + Ep/Eg), 
Ep is the Fermi energy, and we have used the notation 

En,± = {n + l/2)nw, ± {ll2)PogH. 
The density of the holes is given by^ 
{2mHY'^eH ■ 



N'\Ep,H) = 



J2 VEo -Ep-E, 



, (9) 

The position of the Fermi energy Ep{H) is determined 
by the condition of charge neutrality: 

3 

Z = ^N^'{Ep,H)-N'\Ep,H) = 0. (10) 

i=l 

Below we shall apply the formulas of this section to the 
case when the magnetic field is directed along the trigonal 
axis of the crystal, i.e., when h = (0, 0, 1). 



III. NERNST COEFFICIENT 

Let the z axis of the coordinate system be directed 
along the trigonal axis of the crystal, while the x and y 
axes lie in the basal plane of bismuth. In the presence of a 
weak electric field E and of a small temperature gradient 
VT, the charge-current density j is determined by 



ji = o-ikEl - aik'^kT, 



(11) 



where the subscripts i, k stand for the spatial coordinates 
x,y,z; E* = E + |V/i; ^ is the local chemical potential in 
the crystal; and are the conductivity and Peltier 
(thermoelectric) tensors, respectively. In the absence of 
an external magnetic field the tensors (Tik and are 
diagonal, and one has 



''VV 



(12) 



If the external magnetic field H is applied along the z 
axis, the off-diagonal terms axy, (Jyx and axy^ ctyx become 
different from zero, and they satisfy the relationships: 



^xy 



'yx. 



ixy 



J^yx- 



(13) 



In this situation the temperature gradient V^jT will gen- 
erate a transverse electric fie 
ly is defined by the relation 



crate a transverse electric field Ey. The Nernst coefficient 



'WxT' 

and it can be found from the formula 



Jxy 



vH 



■' xx^xy 
(j2 



(14) 



(15) 



xy 



3 



which immediately follows from Eqs. pT|) - (IT4l) . 

Since the oscillations of the Nernst coefficient can be 
observed only at sufficiently high magnetic fields, we as- 
sume below that the condition 



UJcT > 1, 



(16) 



is fulfilled for both the electrons and the holes in bismuth 
where r is the scattering time of an electron or a hole on 
impurities. In clean crystals of bismuth condition (jl6p 
can be satisfied even for moderate magnetic fields. For 
example, using the parameters presented in Table I and 
the experimental data of Ref. we estimate that at 
H = 1 kOe directed along the trigonal axis the values 
of ljJcT for the electrons and holes are of the order of 10 
for the crystal used by K. Behnia et alJ^ (here we have 
assumed that the scattering lengths for the holes and the 
electrons are approximately the same). 

Under condition (1161) . the diagonal component a^x of 
the Peltier tensor is small as compared with axyi^ 



1 



< 1. 



(17) 



Since bismuth has equal numbers of electrons and holes, 
the ratio (Jxy/o'xx in this semimetal is also suppressed 
as compared to the case of a metal with a single type 
of charge carrier™ (e.g., in Ref. [l^ axx was approxi- 
mately ten times larger than axy)- Then, we obtain from 
Eq. ^ that 



(18) 



'xy 



where pik — cr-f, is the resistivity tensor. As the experi- 



mental data show 



is a smooth and almost mono- 



tonic function of the magnetic field in that region of the 
T — H plane in which the oscillations of the Nernst coef- 
ficient were observed f^i^ Hence, the peaks in the Nernst 
coefficient mainly result from the peaks in axy, and be- 
low we shall analyze only this axy, taking Pxx{H) from 
the experiment.— 

The off-diagonal term of the Peltier tensor, axy, con- 



carriers in bismuth, 



dZ{E, H) 
'dE 



E 



dNf{E,H) 
dE 



dN^{E,H) 



dE 



(21) 



According to formulas ([T]), ([6]), ([U, and ([9]), the deriva- 
tive dZ/dE diverges when the energy E crosses the Lan- 
dau levels. The integral over E with the hat-like function 
{fp ^nfp + (1 — fp) ln(l — fp)} of characteristic width 
T smooths the singularities of dZ/dE, but at low tem- 
peratures, fiwc ^ T, the quantity aiy^ shows sharp peaks 
when the Fermi energy approaches a Landau level. These 
simple considerations qualitatively explain the origin of 
the peaks in the Nernst coefficient of bismuth. 

The phonon-drag part a^'xy^ originates from the charge 
carriers that are dragged by a phonon flow in the crys- 
tal due to the electron-phonon interaction. In the 
strong magnetic fields this part is described by a suffi- 
ciently complicated formula which schematically has the 



a^P^'^ oc — 



(22) 



where Vp^ is the probability of phonon scattering by the 
charge carriers; i/o is the probability of phonon scatter- 
ing by other phonons and by the boundaries of the sam- 
ple; (. . .) means some averaging over the wave vectors of 
the phonons existing at the temperature T. It is essen- 
tial here that the probability is proportional to the 
density of states dZ/dE and oscillates with changing H . 
At very low temperatures the phonon scattering by the 
boundaries of the sample prevails over other phonon scat- 
tering processes,— i.e., one has va ^ Vpe. In this case the 

expression for a^iy"^ is proportional to dZ /dE and is sim- 
ilar to Eq. (I^ni) . In order to estimate the magnetic-field 
dependence of the phonon-drag part in this situation, we 
shall use the simplest approximation 



a 



iph) 

xy 



{H,T) 



r{T)a(^{H,T), 



(23) 



sists of the purely electronic and the phonon-drag parts, and axy at sufficiently low magnetic fields when the os- 



where the factor r(T) can be found as the ratio of ai'j/''' 

(e) 



ax 



xy 



xy ■ 



(19) 



The purely electronic part a^cy is caused by a nonequilib- 
rium distribution of the electrons and holes over the crys- 
tal when VxT ^ 0. In strong magnetic fields, uJcT ^ 1, 
a\cy is dissipationless, i.e. it does not depend on the re- 
laxation mechanisms of this distribution>i^ii^ 



cillations of axy become small. With increasing T the 
probability Vpe may become of the order of or even 
exceed it. According to Ref. lo^, this seems to occur in 
bismuth at T ~ 3 — 4 K. Then, Vpp/ {vpe + ^^o) Ij and 
the oscillations of dZ/dE with changing H scarcely re- 
veal themselves in a 



(ph) 



cS 



c 

H 



+00 f)y 

dEg^ifp Infp + (1 - fp) HI - fp)}, 



(20) 

where S is the entropy of the charge carriers per a unit 
volume; fp is the Fermi function with the chemical po- 
tential fj,, and dZ/ dE is the density of states of the charge 



Effect of impurities on the electron and hole 
densities of states in bismuth 



Although in the limit of strong magnetic fields, oJcT ^ 
1 , formula ([20)1 does not depend on t explicitly, the scat- 
tering of the charge carriers on impurities renormalizes 
the density of states of the electrons and holes in bismuth. 
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causing a smearing of the peaks in axy We calculate the 
renormalized dZ/dE for bismuth using the approach of 
Ref . M- 

Within this approach Green functions Gi of charge car- 
riers in the magnetic field H are considered that are aver- 
aged over the positions of impurities in the crystal. The 
index I means the set of quantum numbers n, pz, Py, 
Sj, = ±1 that characterize eigenstates of a charged parti- 
cle in a magnetic field in the absence of impurities. The 
potential of an impurity, U{r), is assumed to be of the 
form U{r) ~ /<5(r) where 6{r) is the delta function, and 
the constant factor / determines the scattering ampli- 
tude of a charged particle on an impurity at e = and 
H ~ 0. The density of states of the charge carriers per 
unit interval of their energy e, dN{e)/de, is expressed in 
terms of the Green functions as 



^ = ilniG(e), 
oe TT 



(24) 



where G{e) = Gi{e). The function G{£) can be found 
from the self-consistent equatioi*^^ 



£l 



(25) 



where £i = s{n,pz, Sz) is the spectrum of a charge carrier 
in the magnetic field in the absence of impurities, and 
Ci is the impurity concentration. The quantity / can 
be expressed in terms of the scattering time r of the 
electrons (the holes) at _ff = when the energy e of these 
particles is equal to the chemical potential fio of bismuth 
in zero magnetic field: 



iPG\H~ 



^27- 



(26) 



Below we calculate G under the simplification that the 
scattering of the holes and of the electrons of the three 
different ellipsoids of bismuth does not mix these parti- 
cles, and so one has electron Tg and hole Th scattering 
times. 

The integration over pz and py in Eq. (j25p can be done 
explicitly for the electrons and holes that are described 
by the spectra ^ and respectively. Then, in the 
case of the holes we obtain the following equation in G: 



G-Go 



where 



1/2 TT 

. ■m,z en 



9n,±{x) 



Y.{9nM'G)-9n,±{^)}, (27) 



(28) 



and Go{e) is the function G{e) in the absence of impuri- 
ties. The function Go is imaginary and is given by 



ImGo 



de 



(29) 




FIG. 1. The density of states for the electrons in bismuth. 
The solid line has been calculated with Eqs. (|24|| - (|3Hl and 
h/Te = 0.001 meV. The points show the analytical result ob- 
tained with the substitution (|32p at h/r^ =0.1 meV. 
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FIG. 2. The density of states for the holes in bismuth. The 
solid line has been calculated with Eqs. (I24|| - (|3ip and fi/r^ = 
0.5-10"'^ meV. The points show the analytical result obtained 
with the substitution ^ at h/rl^ = 0.007 meV. 



where ONq j de is the density of states of the holes in the 
absence of impurities, i.e, dN^ /de directly follows from 
Eq. ([9]). The functions G for all the electron groups are 
identical, and the equation in G has the same form as 
Eq. (HZl) but with 



9n,±{x) 



and 



Eg -f 2(e + x) 



eI'^ [{e + xY -f Eg{e + x- £;„,±)] 



1/2 



ImGo 



de 



, (30) 



(31) 



where dN^/de is the density of states of one of the elec- 
tron ellipsoids in absence of the impurities. The renor- 
malized densities of states of the electrons and of the 
holes in bismuth calculated with equations ([24])- (|3T|) are 
shown in Figs. 1 and 2. 

Interestingly, the renormalized densities of states can 
be approximately obtained by the following simple sub- 
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stitution in formulas ([8]), ([9]): 



+ r2 



(32) 



and by the subsequent differentiation of these formulas 
over Ep. Here T — H/t' , u — EQ — Ep — En.± for the holes 
and u = Ep — for the electrons. It turns out that 

at an appropriate choice of the effective scattering times 
r', the densities of state thus obtained very well describe 
the densities of states calculated with Eqs. (|24l)- ([3T|) . 
see Figs. 1 and 2. The effective scattering time r' weakly 
depends on the number n of the Landau level nearest to 
Ep, and so one and the same value of r' can describe 
the oscillations of the density of states in a rather wide 
interval of the magnetic fields, Figs. 1 and 2. Note that 
r' found for strong magnetic fields is essentially larger 
than T defined in the limit — >■ 0. 

The substitution ([32]) can be understood from the fol- 
lowing considerations: Let us average the density of 
states of charge carriers in the absence of impurities, 
dNo{E)/dE, with a hat-like function characterized by 
some half-width F = H/t': 



dNje) 



I 



dE 



ih/r') 



dNo{E) 



{E - e)2 + {h/r'y dE 



(33) 



In the case of the holes in bismuth (and for any quadratic 
spectrum) this integral is calculated analytically, and the 
density of states obtained coincides with dN{Ep)/dEp 
calculated with the substitution (|32|) . 

Using the substitution, one can easily understand the 
result that follows from the calculations with Eqs. 
pip : In the presence of impurities, the positions of the 
peaks in dN{Ep)/dEp generally do not coincide with 
the magnetic fields at which the Landau levels calculated 
at zero impurity concentration, Ci = 0, cross the Fermi 
energy. Indeed, for the holes in bismuth, we find from 
formula ([9]) with substitution (|32]) that 



dEp 



E 



F2 



Un,± 



F2 



(34) 

where Un,± — Eq — Ep — E„^±. Now each term of this 
sum reaches its maximum at Un,± = On the other 

hand, in the absence of impurities the singularities in 
dNo{Ep) /dEp occur at Un.± = 0. In other words, if one 
takes into account the impurities, the positions of the 
peaks in the density of states shift. 

IV. RESULTS 

In Figs. [SHl] we show the magnetic-field dependences 
of the quantity S^y = vH calculated with Eqs. (IT51) - ([31]) 
at different temperatures T, using the data on the resis- 
tivity Pxx{H) measured in Ref. !l4'. For comparison, in 
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FIG. 3. The Nernst signal S^y = i^H as a function of the 
magnetic field, H, ai T = 0.18 K. The points show the exper- 
imental data,— and the solid line is the result of the calculation 
with Eqs. (UHl) - dSTJ. Here r(0.18) = 0.85, H/th = 0.005meV, 
and h/Te — 0.005meV. The arrows indicate positions of the 
electron peaks visible in the solid line. 




3 4 5 6 

H (10 kOe) 



8 10 



FIG. 4. As Fig. E] but here T = 0.49 K and r(0.49) = 2.5. 



these figures we also present the appropriate experimen- 
tal data^ for Sxy[H). In the calculation of the curves 
Sxy{H) we take into account the i?-dependence of the 
chemical potential /i which is determined by Eq. (jlOp . 
Since the renormalization of the electron and hole den- 
sities of states by impurities influences this dependence, 
we have calculated and in equation using 
formulas ([8]), ^ and substitution (|32|) . At tempera- 
tures T = 0.18 K, 0.49 K and 1.04 K the factor r in 
Eq. (|23p has been chosen in such a way that the calcu- 
lated and the experimental curves coincide at a point near 
the lowest magnetic field H ^ 10 kOe, i.e. at the point 
where the oscillations of Sxy(H) are small. At T = 4.3 K 
(Fig. [5]), in accordance with the considerations of Sec. Ill, 
we take si^^\H) as a constant. The value of this con- 
stant, si^^^ — 0.21 mV/K, has been found from the con- 
dition that the calculated Sxy at = 10 kOc fits the 



appropriate experimental value. 

An inspection of Figs. ISHHl shows that the theoretical 
curves sufficiently well describe the experimental data 
especially in the low-field part of the figures where the 
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FIG. 5. As Fig. Elbut here T = 1.04 K and r(1.04) = 4.5. 



electron peaks essentially do not manifest themselves. 
Note also that the values Th — ^ 1.2 • 10^^*^ sec used 
in the construction of these figures are of the same or- 
der of magnitude as those estimated in Ref. [13 from the 
data of transport measurements (r/j w 8.5 • 10~^° sec and 
Te ~ 4.4- 10~^° sec). Interestingly, although is equal to 
Th here, we find a relatively small amplitude for the elec- 
tron peaks since the hole density of states is essentially 
larger than the electron contribution to dZ{Ep)/dEp. 
That is why only a few electron peaks are seen in the 
experiment,— and the value of Tg has a weak effect on the 
calculated Sxy The visible difference in the positions of 
the electron peaks in the theoretical and experimental 
curves seems to be caused by an insufficient accuracy of 
the spectrum of Smith et aZ.— used here. As to the dis- 
crepancy between the hole peaks in the theoretical and 
experimental curves in the high-field part of the figures, 
it could be due to the two reasons. First, our simplest 
approximation for the phonon-drag part a^^^ , Eq. ([23| , 
and for sif{T = 4.3K) [i.e., sif^T = 4.3K) = 0.21 
mV/K] is probably impaired as magnetic field H devi- 
ates farther and farther away from the point 10 kOe at 
which the factors r and the value of s'^y^\T = 4.3K) 
were found. Second, equation (1^ obtained in Ref. [13, 
strictly speaking, is valid in the semiclassical limit (when 
the Landau- level number n ^ 1), while we use this equa- 
tion in the whole interval of the magnetic fields. In other 
words, in the ultraquantum limit (n ^ 1) a special anal- 
ysis of Sxy is required. 

In the inset of Fig. IH] we compare the values of Sxy 
calculated at iJ = 1 kOe for T = 0.18 K, 0.49 K, 1.04 
K, and 4.3 K with the experimental temperature depen- 
dence of Sxy measured in Ref. [l3 at the same magnetic 
field. These two functions Sxy{T) are seen to agree rea- 
sonably well if it is remembered that they correspond to 
different crystals of bismuth [the phonon-drag part 
depends on the dimensions of the sample at low temper- 
atures when the phonon mean- free path I cx {vpe + vq)^^ 
is set by the sample size] . 

The values of the factor r obtained (see the captions 
to Figs. [3HS]) demonstrate that the relative contribution 
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FIG. 6. As Fig. E] but here T = 4.3 K and instead of formula 
p3)l we use S^Iy^[H) = 0.21 mV/K. The inset shows the 
function Sxy{T) at H = 1 kOe (the solid line represents the 
experimental data,— while the points are the values calculated 
here) . 

of the phonon-drag part 5^^''' to Sxy decreases with de- 
creasing T. Our estimates show that at T < 0.05 K 
the phonon-drag part will practically die out, and one 
will be able to neglect S'i^'*-' completely in the calcula- 
tions and hence to compare directly the electronic part 
S'iy' (H) with the experimental data. Therefore, the mea- 
surements of the Nernst signal at such low temperatures 
could provide the possibility to compare the theory and 
the experiment more accurately. 



V. CONCLUSIONS 

We calculate the purely electronic part SjfJ (H) of the 
Nernst signal in bismuth at low temperatures for the case 
when the magnetic field is directed along the trigonal axis 
of the crystal. In this calculation we take into account the 
renormalization of the densities of states of the electrons 
and holes in bismuth due to scattering of these charge 
carriers on impurities. Without this renormalization, the 
amplitude of the oscillations of S'>xy{H) at low temper- 
atures would be too large as compared with the ampli- 
tude observed in the experiment.^ As to the phonon-drag 
part s'^^\h) of the Nernst signal, this part has been 
taken into account within the simplest approximation. 
The results of the calculations of Sxy = sl^y -\- si^'^^ are 
compared with the experimental dataj^iiS and we find a 
reasonable agreement between the theoretical and exper- 
imental curves in the region of the usual peaks of the 
Nernst coefficient {H < 120 kOe). We also explain the 
lesser importance of the electrons than the holes in the 
production of the peaks and point out that measurements 
of the Nernst signal at lower temperatures would be use- 
ful for a more accurate comparison of the theory and 
the experiment. The approach of this paper, in princi- 
ple, can also be used in analyzing the Nernst-coefficient 
oscillations recently observed in graphitci^ 
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